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ABSTRACT 

We present photometry of two transits of the giant planet WASP-4b with 
a photometric precision of 400-800 parts per milhon and a time samphng of 
25-40 s. The two midtransit times are determined to within 6 s. Together with 
previously published times, the data are consistent with a constant orbital period, 
giving no compelling evidence for period variations that would be produced by a 
satellite or additional planets. Analysis of the new photometry, in combination 
with stellar-evolutionary modeling, gives a planetary mass and radius of 1.237 ± 
0.064 Mjup and 1.365 ±0.021 -Rjup. The planet is 15% larger than expected based 
on previously published models of solar-composition giant planets. With data of 
the quality presented here, the detection of transits of a "super-Earth" of radius 
1.75 i?0 would have been possible. 

Subject headings: planetary systems — stars: individual (WASP-4=USNO-B1.0 0479- 
0948995) 



1. Introduction 

Wilson et al. (2008) recently reported the discovery of WASP-4b, a giant planet that 
orbits and transits a G7V star with a period of 1.34 days. This discovery is notable because 
the planet has an unusually large radius and short orbital period, and the star is one of the 
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brightest transit hosts that is known in the Southern sky {V — 12.5). The large radius seems 
to place the planet among the "bloated" planets for which there is no clear explanation (see, 
e.g., Burrows et al. 2007, Guillot 2008). The short period raises the possibility of observing 
tidal decay (Rasio et al. 1996, Sasselov 2003) and makes WASP-4b an attractive target for 
observations of occultations (secondary eclipses) that would lead to detections of the reflected 
light and thermal emission from the planet's atmosphere. 

The host star's brightness and southern declination are important for a practical rea- 
son: they allow the large-aperture telescopes of the southern hemisphere to be used ad- 
vantagcoTisly. In this paper, we report on observations of a transit of WASP-4b with the 
MagcUan/Baadc 6.5m telescope, with the goal of deriving independent and refined param- 
eters for this interesting system. Previous papers in this series, the Transit Light Curve 
project, have achieved this goal by combining the information from many independent tran- 
sit observations with smaller telescopes (see, e.g., Holman et al. 2006, Winn et al. 2007). 
In principle, with a larger telescope, it should be possible to achieve this goal with fewer 
observations and also to measure precise midtransit times, which can be used to search 
for additional planets via the method of Holman & Murray (2005) and Agol et al. (2005). 
Gillon et al. (2008) recently presented photometry of WASP-4 with one of the 8.2m Very 
Large Telescopes, with the same motivation. 

This paper is organized as follows. § 2 describes the observations and data reduction, § 3 
describes the photometric analysis, § 4 describes the results of stellar-evolutionary modeling, 
and § 5 discusses the newly-measured midtransit times and a refined cphcmcris. Finally, § 6 
discusses the refined measurement of the planetary radius, and considers how small a planet 
we could have detected, given data of the quality presented here. 

2. Observations and Data Reduction 

We observed the transits of UT 2008 Aug 19 and 2008 Oct 09 with the Baade 6.5m 
telescope, one of the two Magellan telescopes at Las Campanas Observatory in Chile. We 
used the Raymond and Beverly Sackler Magellan Instant Camera (MagIC) and its SlTe 
2048 X 2048 pixel CCD detector, with a scale of (^.'069 pixel"^. Ordinarily this detector uses 
four amplifiers, each of which reads a quadrant of the array, giving a total readout time of 
23 s. We used a 2048 x 256 pixel subarray and a single amplifier, giving a readout time of 
10 s. We rotated the field of view to align WASP-4 and a nearby comparison star along the 
long axis of the subarray (parallel to the read register). The comparison star is 36" east and 
71" south of WASP-4. At the start of each night, we verified that the time stamps recorded 
by Magic were in agreement with GPS-based times to within one second. On each night 
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we obtained repeated z-band exposures of WASP-4 and the comparison star for about 5 hr 
bracketing the predicted transit time. Autoguiding kept the image registration constant to 
within 10 pixels over the course of the night. 

During the 2008 Aug 19 observations, WASP-4 rose from an airmass of 1.09 to 1.03, 
and then set to an airmass of 1.23. At first, we used an exposure time of 30 s. Shortly 
after midtransit the seeing improved abruptly, from a full-width at half-maximum (FWHM) 
of 11 pixels to 7 pixels. As a result of the higher rate of detected photons per pixel, some 
images were spoiled due to nonlinearity and saturation. For the rest of the night we used an 
exposure time of texp = 15 s. 

During the 2008 Oct 09 observations, WASP-4 rose from an airmass of 1.10 to 1.03, 
and then set to an airmass of 1.21. Having learned our lesson in August, we deliberately 
defocused the telescope so that the image width was dominated by the effects of the telescope 
aberration rather than natural seeing. This was done by moving the secondary mirror by 
a constant displacement relative to the in-focus position determined by the guider probe. 
The exposure time was texp = 30 s. The stellar images were "donuts" with a diameter of 
approximately 25 pixels. 

We used standard IRAF procedures for overscan correction, trimming, bias subtraction, 
and flat-field division. The bias frame was calculated from the median of 60 zero-second 
exposures, and the fiat field for each night was calculated from the median of 60 z-band 
exposures of a dome fiat screen. We performed aperture photometry of WASP-4 and the 
comparison star and divided the flux of WASP-4 by the flux of the comparison star. We 
experimented with different aperture sizes and sky regions, aiming to minimize the variations 
in the out-of-transit (OOT) portion of the differential light curve. Best results were obtained 
with an aperture radius of 38 pixels for the 2008 Aug 19 observations and 35 pixels for the 
2008 Oct 09 observations. 

A few images were also obtained with Johnson-Cousins VRI filters to measure the 
difference in color between WASP-4 and the comparison star. The instrumental magnitude 
differences (target minus WASP-4) were AV = 0.434, AR = 0.469, AI = 0.478, and Az = 
0.472. Evidently the two stars are similar in color, with the comparison star being slightly 
bluer [A(V^ — I) = —0.044]. As described in § 3, the 2;-band time series was corrected for 
differential extinction between the target and comparison star by fitting a linear function of 
airmass to the magnitude difference. The results were consistent with the expectation that 
the bluer comparison star suffers from greater extinction per unit airmass. The extinction- 
corrected data are given in Table [H and plotted in Fig. [H along with the best-fitting model. 

Due to the abrupt seeing change on 2008 Aug 19 and the associated change in exposure 
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Fig. 1. — Top. — Relative z-band photometry of WASP-4 based on observations with the 

Magellan (Baade) 6.5m telescope. Middle. — Composite light curve. The solid line shows the 
best-fitting model. Bottom. — Residuals between the data and the best-fitting model. 
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level, we consider the 15 s and 30 s exposures as two separate time series (TSl and TS2). 
Thus, together with the 2008 Oct 09 time series (TS3), there were 3 time series to be 
analyzed. The TSl pre-ingress data has a median time sampling of 41 s and a standard 
deviation of 478 parts per million (ppm). The TS2 post-egress data has a median time 
sampling of 26 s and a standard deviation of 691 ppm. These noise levels are about 13% 
larger than the calculated noise due to photon-counting statistics, read noise, sky noise, and 
scintillation noise (using the approximate formulas of Reiger 1963 and Young 1967). The 
ratio of the observed TSl noise to the observed TS2 noise is 0.69, which is nearly equal to 
(^exp,2/^exp,i)~^'^^ = 0.71. The near-equality is evidence that the dominant noise source is a 
combination of photon noise and scintillation noise, both of which vary as Uxp^- In the final 
version of TS3, which has a median time sampling of 41 s, the pre-ingress data has an rms of 
475 ppm and the post-egress data has an rms of 488 ppm. These values are about 17% above 
the calculated noise level. In neither light curve did we detect signficant correlation between 
the noise and the pixel coordinates, FWHM, or shape parameters of the stellar images. 



3. Photometric Analysis 

We fitted a transit light curve model to the data, based on the analytic formulae of 
Mandel & Agol (2002). The set of model parameters included the planet-to-star radius ratio 
{Rp/R-i,), the ratio of the stellar radius to the orbital semimajor axis (R^^/a), the orbital 
inclination {{), and the midtransit time (Tg). We also fitted for three parameters (Amo, kz 
and kt) specifying a correction for systematic errors, 

Amcor = Amobs + Amo + Kz + ktt, (1) 

where z is the airmass, t is the time, Amobs is the observed magnitude difference between 
the target and comparison star, and Amcor is the corrected magnitude difference that is 
compared to the idealized transit model. The kz term specifies the differential airmass 
extinction correction that was mentioned in § 2. The kt term was not strictly necessary 
to fit the data (it was found to be consistent with zero), but we included it in order to 
derive conservative error estimates for the transit times, as the parameter kf is covariant 
with the transit time. We allowed Amo, kt, and kz to be specific to each time series, with 
the exception that TSl and TS2 had a common value of kz because those data were obtained 
on the same night. 

The final two model parameters were the coefficients {ui and U2) of a quadratic limb- 
darkening law, 

^ = l-«i(l-/x)-«2(l -//)', (2) 
-'1 
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where ji is the cosine of the angle between the hne of sight and the normal to the stellar 
photosphere, and is the specific intensity. We allowed Ui and U2 to vary freely subject 
only to the conditions Ui + U2 < 1, Ui + U2 > 0, and Ui > 0. It proved advantageous to 
perform the fit using the linear combinations 

vi — ui + 2.33 1*2, V2 — ui — 2.33it2, (3) 

because vi and V2 have nearly uncorrelated errors (for further discussion, see Pal 2008). 

We assumed the orbit to be circular because the RV data of Wilson et al. (2008) are 
consistent with a circular orbit, and because the expected timescale for orbital circularization 
at present. 




(Rasio et al. 1996) is 0.3 Myr (Qp/lO^) for WASP-4b, which is much shorter than the 
estimated main-sequence age of the star. In this expression, Qp is the the tidal dissipation 
parameter (see, e.g., Goldreich & Soter 1966). This order-of-magnitude argument suggests 
that assuming a circular orbit is reasonable, although the value of Qp for irradiated giant 
planets is highly uncertain, and the expected timescale is highly approximate because it 
ignores the coupled evolution of the orbital distance and eccentricity (Jackson, Greenberg, 
& Barnes 2008). 

The fitting statistic was 

N 
1=1 

where N is the number of flux measurements (photometric data points), /j(obs) is the ith 
measurement, /^(calc) is the calculated flux given a particular choice of model parameters, 
and (7/,i is the uncertainty in the ith measured flux. We determined appropriate values of 
af^i as follows. First, we multiplied the calculated errors in each time series by a constant 
chosen to give a minimum value of Xpf-^dof — 1- The constants were 1.20, 1.12, and 1.21 
for TSl, TS2, and TS3 respectively. Next, we assessed the time-correlated noise (also cafled 
"red noise") by examining the autocorrelation function, the power spectrum, and a plot of 
the Allan (1966) deviation of the residuals. Wc also used the method described by Winn 
et al. (2008), in which the ratio /3 is computed between the standard deviation of time- 
averaged residuals, and the standard deviation one would expect assuming white noise. 
For TSl and TS2 we found no evidence for signiflcant correlations. For TS3 we found 
structure in the autocorrelation function on a time scale of 15-20 min, the approximate 
ingress or egress duration, giving P — 1.52. One naturally suspects that the correlated noise 
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represents measurement error, although it is also possible that the noise is astrophysical, 
arising from starspots or other stellar inhomogeneities. In support of an astrophysical origin, 
we find no evidence for correlated noise {(5 = 1) when considering only the out-of-transit 
data. Nevertheless we cannot draw a firm conclusion; instead we attempt to account for 
the correlations by multiplying the error bars of TS 3 by an additional factor of f3 = 1.52. 
Thus, for TSl, TS2, and TS3, the final values of af^i were equal to the calculated error bars 
multiplied by 1.20, 1.12, and 1.21 x 1.52 = 1.84, respectively. The error bars given in Table 1 
are the final values of cr/^j that were used in the fitting process. 

To determine the allowed ranges of each parameter, we used a Markov Chain Monte 
Carlo (MCMC) technique, with the Gibbs sampler and the Metropolis-Hastings algorithm, 
to estimate the a posteriori joint probability distribution of all the model parameters. This 
algorithm creates a sequence of points (a "chain") in parameter space by iterating a jump 
function, which in our case was the addition of a Gaussian random deviate to a randomly- 
selected single parameter. After this operation, if the new point has a lower Xf than the 
previous point, the "jump" is executed: the new point is added to the chain. If not, then 
the jump is executed with probability exp {-Axf/2). When the jump is not executed, the 
current point is repeated in the chain. The sizes of the random deviates are adjusted so 
that ~40% of jumps are executed. After creating 10 chains of 500,000 links to check for 
mutual convergence, and trimming the first 20% of the links to eliminate artifacts of the 
initial conditions, the density of the chains' points in parameter space was taken to be the 
joint a posteriori probability distribution of the parameter values. Probability distributions 
for individual parameters are created by marginalizing over all other parameters. 

The results are given in Table [21 which gives the median of each distribution, along with 
the 68.3% lower and upper confidence limits (defined by the 15.85% and 84.15% levels of 
the cumulative distribution). The entries designated A are those that follow directly from 
the photometric analysis. The entries designated B are those that are drawn from Gillon et 
al. (2008) and are repeated here for convenience. The entries designated C are based on a 
synthesis of our modeling results and theoretical models of stellar evolution, as discussed in 
§ 4. The entries designated D are the parameters of a refined transit ephemeris based on the 
two newly- measured transit times as well as some other available timing data (see § 5). 

As a consistency check we tried fitting the 2008 Aug 19 data and the 2008 Oct 09 data 
separately. We found that the results for the parameters Rp/R^,, i?*/a, i, ui, and U2 were all 
in agreement within la, suggesting that our error estimates are reasonable. 
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4. Theoretical isochrone fitting 

The combination of transit photometry and the spectroscopic orbit (radial- velocity vari- 
ation) of the host star do not uniquely determine the masses and radii of the planet and the 
star. There remain fitting degeneracies Mp oc M^^^ and Rp oc (x M^^ (see, e.g., Winn 
2008). We broke these degeneracies by requiring consistency between the observed properties 
of the star, the stellar mean density that can be derived from the photometric parameter 
a/R^ (Seager & Mallen-Ornelas 2003, Sozzetti et al. 2007), and theoretical models of stellar 
evolution. The inputs were T^s = 5500 ± 100 K and [Fe/H] = -0.03 ± 0.09 from Gillon et 
al. (2008), the stellar mean density = 1.694lQQ3y g cm~^ derived from the results for the 
a/R^ parameter, and the Yonsei-Yale (Y^) stellar evolution models by Yi et al. (2001) and 
Demarque et al. (2004). We computed isochrones for the allowed range of metallicities, and 
for stellar ages ranging from 0.1 to 14 Gyr. For each stellar property (mass, radius, and 
age), we took a weighted average of the points on each isochrone, in which the weights were 
proportional to exp(— x^/2) with 
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(6) 



Here, the A quantities denote the deviations between the observed and calculated values 
at each point. The asymmetric error bar in p^ was taken into account by using different 
values of depending on the sign of the deviation. The weights were further multiplied by 
a factor taking into account the number density of stars along each isochrone, assuming a 
Salpeter mass function. This procedure is essentially the same as that employed by Torres 
et al. (2008). The only difference is that we calculated the 68.3% uncertainties by assuming 
that the errors in Teg, p*, and [Fe/H] obey a Gaussian distribution, while Torres et al. (2008) 
took the distribution to be uniform within the quoted la limits. 

Through this analysis, we found = 0.925±0.040 Mq and R^ = 0.912±0.013 Rq. The 
stellar age was poorly constrained, with a formally allowed range of 6.5±2.3 Gyr and a nearly 
uniform distribution. The corresponding planetary mass and radius were obtained by merg- 
ing the results for the stellar properties with the parameters determined in our photometric 
analysis, and with the stellar radial-velocity semiamplitude = 0.24±0.01 km s~^ measured 
by Wilson et al. (2008). The results are Mp = 1.237±0.064 Mj^p and Rp = 1.365±0.021 i?jup. 
Table [2] gives these results, along with the values for some other interesting parameters that 
can be derived from the preceding results. As a consistency check, we computed the implied 
stellar surface gravity and its uncertainty based on our analysis, finding logs' = 4.481 ±0.008 
where g is in cm s~^. This agrees with the spectroscopic determination of surface gravity, 
log (7 = 4.3 ± 0.2, made by Wilson et al. (2008) based on an analysis of the widths of 
pressure- sensitive lines in the optical spectrum. 
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It is important to keep in mind that the quoted error bars for the parameters designated 
C in Table [2] are based on the measurement errors only, and assume that the any systematic 
errors in the isochrones are negligible. As a limited test for the presence of such errors, 
Torres et al. (2008) tried analyzing transit data for 9 systems using isochrones computed by 
three different groups: the isochrones used here as well as those of Girardi et al. (2000) 
and Baraffe et al. (1998). They found the differences to be smaller than the error bars, 
especially for stars similar to the Sun such as WASP-4. Similar results were found by 
Southworth (2009). Nevertheless, the true systematic errors in the Y^ isochrones are not 
known, and we have not attempted to quantify them here, although it seems plausible that 
the masses and radii are subject to an additional error of a few percent. 



5. Transit Times and a Refined Ephemeris 

Precise measurements of transit times are important because the gravitational perturba- 
tions from other bodies in the system (such as a satellite or additional planet) could produce 
detectable variations in the orbital period. Based on the MCMC analysis described in the 
previous section, the uncertainties in the two transit times are 4.6 s and 4.9 s, making them 
among the most precise such measurements that have been achieved. This is a consequence 
of good photometric precision and fine time sampling, along with the relative insignificance 
of time-correlated noise and the large transit depth. 

As a check on the error bars, we also estimated the midtransit time and its error using the 
residual permutation (RP) method, a type of bootstrap analysis that attempts to account for 
time-correlated errors. Fake data sets are created by subtracting the best-fitting model from 
the data, then adding the residuals back to the model after performing a cyclic permutation 
of their time indices. For each fake data set, Xf is minimized as a function of Tc, Amo, kz 
and kt (the other parameters being held constant, as they are nearly uncorrelated with Tc). 
The distribution of the results for Tc is taken to be the a posteriori probability distribution 
for Tc. For the 2008 Aug 19 data, the RP-based error bar was only 4% larger than the 
MCMC-based error bar; the results were nearly indistinguishable. For the 2008 Oct 09 data, 
the RP-based error was 25% larger than the MCMC-based error. To be conservative, we 
report the larger RP-based errors in Table 1 and used those larger error bars in recomputing 
the ephemeris (see below). 

As a further check on the analysis, we allowed TSl and TS2 to have independent values 
of Tc, Amo, kz, and kt during the fitting process. In other words we fitted all the data 
but did not require that TSl and TS2 agree on the transit time. The result was that the 
difference between Tc(TSl) and Tc(TS2) was 2.3 ±13.4 s. We also tried a similar experiment 
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with the 2008 Aug 19 data, sphtting it into two equal parts that were fitted jointly with TSl 
and TS2 (which in this case were required to agree on the midtransit time so as to provide 
a constraint on the transit duration). The result in this case was AT^ = 4.6 ± 11.9 s. The 
mutual consistency of the results suggests that our error bars are reasonable. 

The transits of 2008 Aug 19 and 2008 Oct 09 were separated by 38 orbital periods. By 
calculating ATc/38 we derive an independent estimate of the orbital period, P = 1. 3382369 ± 
0.0000024 days. The most precise determination previously reported was P = 1.3382324 ± 
0.0000029 days, by Gillon et al. (2008), based on 2 years of data. Our period has comparable 
precision, though it is based on only 2 transits separated by 50 days. The difference between 
the two independent period determinations is 0.39 ± 0.32 s. 

We fitted a linear function of epoch to the two newly-measured midtransit times along 
with the 5 transit times given in Table 2 of Gillon et al. (2008). The fit had x^=7.8 with 5 
degrees of freedom. The chance of finding a value of this large by chance is 17%, using the 
quoted la error bars and assuming the errors obey a Gaussian distribution. We deem this an 
acceptable fit, and conclude that the available data do not provide compelling evidence for 
any departures from a constant period. The refined transit ephemeris is Tc{E) = Tc{0) + EP, 
with 

Te(0) = 2, 454, 697.797562 ±0.000043[BJD], 

P = 1.33823214 ±0.00000071. (7) 

Fig. 2 shows a plot of the differences between the observed and calculated transit times. 



6. Discussion 

Our results for the planetary and stellar properties are in accord with the previous 
analyses by Wilson et al. (2008) and Gillon et al. (2008). In general our error bars are 
comparable in size or smaller than those of Gillon et al. (2008), who observed a transit 
with one of the Very Large Telescopes. This consistency is reassuring, especially since our 
error estimates are more conservative in some respects. We have accounted for uncertainty 
in the limb-darkening coefficients, as well as the slopes of systematic trends with time and 
airmass. Previous investigators assumed that these parameters were known exactly, leading 
to underestimated errors in any covariant parameters. Southworth (2008) demonstrated this 
effect for the limb-darkening coefficients in particular^ In addition, Gillon et al. (2008) 



^Our results for the limb-darkening coefficients are Ui ~ 0.311 ± 0.041 and U2 ~ 0.227 ± 0.089. These are 
not too far from the tabulated values 0.265 and 0.303 given by Claret (2004) for a star with the observed 
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Table 1. Relative Photometry of WASP-4 



UcLl V V-Cllljl llj J Lliidll i.-'dLrt; 


T-? Ti 1 1 V 


T TtI Q 1 TT^""\r 


2454697.708386 


1.00005 


0.00050 


2454697.708861 


1.00018 


0.00050 


2454697.709337 


1.00003 


0.00050 


2454697.709810 


0.99975 


0.00050 


2454697.710285 


0.99930 


0.00050 


2454697.710759 


0.99996 


0.00049 


2454697.711234 


0.99917 


0.00050 


2454697.711708 


1.00067 


0.00050 


2454697.712185 


1.00039 


0.00050 


2454697.712657 


1.00066 


0.00049 


2454697.713132 


1.00062 


0.00050 


2454697.713606 


0.99982 


0.00050 


2454697.714082 


0.99976 


0.00049 



Note. — The time stamps in Column 1 represent the 
Barycentric Julian Date at the time of midexposure. The 
uncertainty in Column 3 is the quantity cj- ^ discussed in § 3. 
We intend for this Table to appear in entirety in the electronic 
version of the journal. An excerpt is shown here to illustrate 
its format. The data are also available from the authors upon 
request. 
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Fig. 2. — Transit timing residuals for WASP-4b. The calculated times (using the ephemeris 
given in Eq. [7]) have been subtracted from the observed times. 



Table 2. System Parameters of WASP-4b 



Parameter Value 68.3% Conf. Limits Comment 



Transit ephemeris: 



Reference epoch [BJD] 


2454697.797562 


±0.000043 


D 


Orbital period [days] 


1.33823214 


±0.00000071 


D 


Transit parameters: 








Midtransit time on 2008 Aug 19 [BJD] 


2454697.797489 


±0.000055 


A 


Midtransit time on 2008 Oct 09 [BJD] 


2454748.650490 


±0.000072 


A 


Planct-to-star radius ratio, Rp/Ri, 


0.15375 


-0.00055, +0.00077 


A 


Orbital inclination, % [deg] 


00. ou 


— u.4!U, ^u.yo 


\ 


Scaled semimajor axis, a/i?* 


5.473 


-0.051, +0.015 


A 


Transit impact parameter, b = acosi/Ri, 


0.143 


-0.098, +0.038 


A 


Transit duration [hr] 


2.1660 


-0.0062, +0.0054 


A 


Transit ingress or egress duration [hr] 


0.2913 


-0.0002, +0.01006 


A 


Linear limb-darkening coefficient, ui 


0.311 


±0.041 


A 


Quadratic liinb-darkening coefficient, U2 


0.227 


±0.089 


A 


Semimajor axis [AIJ] 


0.02340 


±0.00060 


C 


Planet-to-star mass ratio, Mj,/M* 


0.00127 


±0.00012 


C 


Stellar parameters: 








Mass, M* [Mq] 


0.925 


±0.040 


C 


Radius, i?* [R©] 


0.912 


±0.013 


C 


Surface gravity, log 3* [cgs] 


4.4813 


±0.0080 


C 


Mean density, p* [g cm~^] 


1.728 


-0.047, +0.016 


A 


Eflfective temperature, Tgff [K] 


5500 


±100 


B 


Metallicity, [Fe/H] 


-0.03 


±0.09 


B 


Projected rotation rate, tJsini* [km s~^] 


2.0 


±1.0 


B 


Luminosity [Lq] 


0.682 


±0.065 


C 


Absolute V magnitude 


5.30 


±0.13 


C 


Age [Gyr] 


6.5 


±2.3 


C 


Planetary parameters: 








Mp [Mjup] 


1.237 


±0.064 


C 


Rp [R-jup] 


1.365 


±0.021 


C 


Surface gravity, g-p [ui s^'^] 


16.41 


±0.75 


A 


Mtuui (l('iisii\'. fip g cm \ 


O.CiOl 


±0.012 


(' 



Note. — (A) Based on the analysis of the new light curves (see § 3). (B) From Gillon et al. (2008). 
(0) Functions of group A and B parameters, supplemented by theoretical isochrones (see § 4). 
For (C) parameters, the error bars are based only on the measurement errors and do not account for 
any possible error in the isochrones. (D) Based on the two newly-measured transit times, as well 
as all entries in Table 2 of Gillon et al. (2008). 
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found a much stronger influence of correlated noise on the results for the midtransit time. 
In the language of § 3, their RP-based error bar was 300%-400% larger than the MCMC- 
based error bar, as compared to 4-25% in our case. We also question the apphcability of 
the "small-planet" approximation used by Gillon et al. (2008). With our data, using the 
small-planet approximation would have led to an erroneous 2.3a shift in Rp/R^,. 

For the key parameter Rp, Wilson et al. (2008) found 1.4161°;°^^ Rj^p, Gillon et al. (2008) 
found 1.304lg;Q42 -Rjup, and we find 1.365 ±0.021 -Rjup. Here it is important to reiterate that 
the estimate of Rp relied on the isochrones, which probably contribute an additional 
systematic error of a few percent. Despite this, it seems clear that WASP-4b is a member 
of the class of "bloated" planets, by which we mean planets that are larger than predicted 
according to theoretical models of solar-composition gas giant planets, even after accounting 
for the intense irradiation from the parent star. 

Bodenheimer et al. (2003) predicted the radii of giant planets as a function of the age, 
mass, and equilibrium temperature, defined as 

1/4 

= (1677 K) (1 - Ay/\ (8) 

where A is the Bond albedo, is the stellar luminosity, a is the Stefan-Boltzmann constant, 
and a is the semimajor axis. In the latter equality we have evaluated Teq for WASP-4b using 
the results given in Table [2l As long as the albedo is not very close to unity, Bodenheimer 
et al. (2003) predict a planetary radius of 1.13 -Rjup for a solar composition at 4.5 Gyr. This 
is smaller than the observed radius by llcr. 

Fortney et al. (2007) presented theoretical models parameterized by mass, age, and an 
effective orbital distance, defined as the distance from the Sun where a hypothetical planet 
would receive the same flux as the actual planet, 

/r \-i/2 

ae = a Ij^j = 0.0281 AU, (9) 

where again we have evaluated the expression as appropriate for WASP-4b. For a solar 
composition at 4.5 Gyr, Fortney et al. (2007) predict -Rp = 1.16 -Rjup, about lOcr smaller 
than the observed value. Neither of these sets of models takes into account the "transit 



effective temperature and surface gravity. However, with data this precise, one should not fix the limb- 
darkening coefficients at the tabulated values. The tables have unquantified errors and these errors are 
strongly correlated with the other photometric parameters. For example, for HD 209458 the tabulated 
values can be ruled out with >99.9% precision (Southworth 2008). 



(1 - A)L, 



-14- 



radius effect" of Burrows et al. (2003) — the underprediction of transit radii by modelers who 
use an inappropriately high pressure to define the radius of their models — but this effect is 
expected to be only a few percent according to Fortney et al. (2007). 

One might suppose that the system is young, and has therefore not had as much time 
to contract from its presumably hot and distended initial state. For an age of 300 Myr, 
Fortney et al. (2007) predict a radius of 1.30 -Rjup, only 3(7 smaller than the observed value 
(and perhaps consistent, given the systematic errors in the isochrones and the transit radius 
effect). However, such a young age is disfavored by our isochrone analysis, and Wilson et 
al. (2008) have also argued that the star is older than about 2 Gyr based on a nondetection 
of Li 1 (6708 A) in the spectrum. 

Another possible resolution is to invoke an extra source of internal heat (of unspecified 
origin) within the planet. Liu et al. (2008) have provided useful fitting formulas to gauge how 
much extra power is required, for a solar-composition planet of a given mass and effective 
orbital distance. They are parameterized by the variable x, defined as the ratio of the extra 
power to the incident power from the host star's radiation. By interpolating Table 3 of Liu et 
al. (2008) we find that to match the observed radius of WASP-4b, one needs an "anomalous" 
heating rate of ~0.15% of the incident power from the host star, or ~8 x 10^^ erg s~^. 
According to the same tables, the equilibrium radius is achieved in ~30 Myr. 

Of course the preceding calculation does not solve the problem of the bloated planet. It 
merely quantifies the power requirement for the unknown heat source. The calculations of 
Liu et al. (2008) were in the context of tidal heating due to orbital circularization. According 
to their Eq. (A3), an orbital eccentricity of 0.002 (Qp/10^) would suffice to produce enough 
power to inflate the planet to the requisite degree. Hence for Qp between 10^ and 10^, the 
required eccentricity for this hypothesis is about 0.002 to 0.02. This is small enough to be 
compatible with the RV data, but larger than expected based on the estimated stellar age 
and the order-of-magnitude timescale for tidal circularization (see § 3), unless there is a third 
body whose gravitational influence is preventing circularization. A very precise measurement 
of the time of occultation of WASP-4b by its host star might allow the orbital eccentricity 
(or more precisely e cos a;, where cu is the argument of pericenter) to be measured well enough 
to test this notion. Of course, many other ways to explain the radii of the bloated planets 
have been given in the literature (Guillot & Showman 2002, Showman & Guillot 2002, 
Bodenheimer et al. 2003, Burrows et al. 2007, Chabrier & Baraffe 2007). This issue remains 
unresolved for the particular case of WASP-4b and for the entire ensemble. 

In addition to refining the estimates of planetary parameters, and seeking additional 
bodies through transit timing, an important application of precise differential photometry 
is checking for transits of planets that have been detected by the Doppler technique but for 
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which the orbital inchnation is unknown. Given the progress of the Doppler surveys in finding 
lower-mass planets, it is also a likely route to the discovery of the first examples of "super- 
Earth" planets that transit bright stars, and that are therefore valuable for understanding the 
sizes, orbits, and atmospheres of terrestrial planets. It is generally thought that detection of 
transits of terrestrial-mass planets must be done from space, to eliminate the contaminating 
effects of Earth's atmosphere and to observe for long intervals without interruptions due to 
bad weather and the day/night cycle. These are clear advantages but they must be weighed 
against the high cost of space missions. 

It is interesting to contemplate how small a planet we would have been able to detect, 
given data of the quality presented here. Our light curve has a signal-to-noise ratio of 
approximately y/TTS/a ~ 920, using the notation of Carter et al. (2008) in which FT is the 
number of data points per transit duration, S is the transit depth, and a the photometric 
precision. Statistically one might expect that we would have achieved a 5a detection of 
a planet that has an area 920/5 = 184 times smaller than WASP-4b, i.e., a planet with 
Rp = 1.1 However, it is not clear whether a statistically unassailable 5a detection 

would really be credible, given the possible presence of systematic trends, correlated errors, 
uncertainty in predicted transit times, and general caution. 

To judge the belicvability of a super- Earth detection, we added to the residuals of Fig. 2 
a model light curve of a transiting super- Earth. All of the parameters of the model were the 
same as the actual WASP-4b parameters except the planetary radius, which was reduced to 
1.75 i?®, the approximate upper limit that is predicted for terrestrial planets. The result, 
shown in Fig. 3, gives a visual impression of what one might realistically expect, based on 
two nights of observations of a 12th magnitude star with a large telescope. 

We fitted this simulated data in nearly the same way that we fitted the WASP-4 data. 
The only difference is that we fixed the values of the limb-darkening parameters and the 
stellar radius, as befits a detection rather than a characterization experiment. We also 
assumed the orbital period was known to within 10% (3.2 hr). The free parameters were 
(Rp/Ri,); i; the two mid-transit times; and the parameters describing systematic trends, 
Amo, kz and kt, for each time series. The results were that (Rp/Ri,) was consistent with 
the input value and had an error of 14%, dominated by the correlation with the poorly- 
constrained orbital inclination. The transit times were recovered within 5 min. The detection 
is visually apparent in the time-binned light curve shown in the lower panel of Fig. 3. One 
would certainly prefer a longer stretch of out-of-transit data, especially if the orbital period 
has a large uncertainty. Nevertheless the impression given is that the detection of transits 
of super-Earths around Sun-hke stars is within reach of ground-based observations. 
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Fig. 3. — Simulated photometric observations of a transiting super-Earth around a 12th 
magnitude G7V star, based on our WASP-4 data. The residuals shown in Fig. 1 were added 
to an idealized transit model in which Rp = 1.75 and all other parameters were the same 
as in the WASP-4 model. The bottom panel shows a time-binned version of the simulated 
data. 



-17- 



We thank Paul Schechter and Adam Burgasser for observing on our behalf on several 
occasions, for this project and related projects, and for helpful discussions on observing 
techniques. Partial support for this work came from NASA Origins grants NNX09AB33G 
(to M.J.H.) and NNG04LG89G (to G.T.). 

REFERENCES 

Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 
Allan, D. W. 1966, Proc. IEEE, 54, 221 

Baraffe, I., Chabrier, G., Allard, P., & Hauschildt, P. H. 1998, A&A, 337, 403 

Bodenheimer, P., Laughlin, G., & Lin, D. N. C. 2003, ApJ, 592, 555 

Burrows, A., Sudarsky, D., & Hubbard, W. B. 2003, ApJ, 594, 545 

Burrows, A., Hubeny, L, Budaj, J., & Hubbard, W. B. 2007, ApJ, 661, 502 

Carter, J. A., Yee, J. C., Eastman, J., Gaudi, B. S., & Winn, J. N. 2008, ApJ, 689, 499 

Chabrier, G., & Baraffe, 1. 2007, ApJ, 661, L81 

Claret, A. 2004, A&A, 428, 1001 

Demarque, P., Woo, J.-H., Kim, Y.-C, & Yi, S. K. 2004, ApJS, 155, 667 
Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 
Gillon, M., et al. 2008, arXiv:0812.1998l 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 

Goldreich, P., & Soter, S. 1966, Icarus, 5, 375 

Guillot, T., & Showman, A. P. 2002, A&A, 385, 156 

Guillot, T. 2008, Physica Scripta Volume T, 130, 014023 

Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 

Holman, M. J., et al. 2006, ApJ, 652, 1715 

Jackson, B., Greenberg, R., & Barnes, R. 2008, ApJ, 678, 1396 



- 18 - 

Liu, X., Burrows, A., & Ibgui, L. 2008, ArXiv e-prints, 805. larXiv:0805. 17331 
Mandel, K., & Agol, E. 2002, ApJ, 580, L171 
Pal, A. 2008, MNRAS, 1043 

Rasio, F. A., Tout, C. A., Lubow, S. H., & Livio, M. 1996, ApJ, 470, 1187 

Reiger, S. H. 1963, AJ, 68, 395 

Sasselov, D. D. 2003, ApJ, 596, 1327 

Seager, S., & Mallen-Ornelas, G. 2003, ApJ, 585, 1038 

Showman, A. R, & Guillot, T. 2002, A&A, 385, 166 

Southworth, J. 2008, MNRAS, 386, 1644 

Southworth, J. 2009, rarXi\^08Tr32771 

Sozzetti, A., Torres, G., Charbonneau, D., Latham, D. W., Holman, M. J., Winn, J. N., 
Laird, J. B., & O'Donovan, F. T. 2007, ApJ, 664, 1190 

Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324 

Wilson, D. M., et al. 2008, ApJ, 675, L113 

Winn, J. N., et al. 2007, AJ, 133, 1828 

Winn, J. N. 2008, ArXiv e-prints, 807, larXiv:0807.4929l 

Winn, J. N., et al. 2008, ApJ, 683, 1076 

Yi, S., Demarque, P., Kim, Y.-C, Lee, Y.-W., Ree, C. H., Lejeune, T., & Barnes, S. 2001, 
ApJS, 136, 417 

Young, A. T. 1967, AJ, 72, 747 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



